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The electron energy and density matrices in molecular systems are convex in respect 
of the number of particles. So that, the chemical descriptors based on their derivatives 
present the hamper of discontinuities for isolated systems and consequently higher order 
derivatives are undefined. The introduction of the interaction between the physical 
domain with an environment induces a coherent structure for the density matrix in 
the grand-canonical formulation suppressing the discontinuities leading to the proper 
definitions of the descriptors. 
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The electron distribution rearrangements in molecular systems under the influen¬ 
ce of external perturbations, internal conversions, conformational changes or reactive 
interactions are closely related to chemical reactivity which is relevant to understand 
the molecular structured’^ The chemical units involved in these processes are interac¬ 
ting moieties, atoms or functional groups which undergo a flux of charge by electron 
exchange between themd’^ This fact induces such domains to possess a non-integer (or 
fractional) number of electrons which may be interpreted as an ensemble average in a 
quantum state of an open system at equilibrium after the rearrangement processd’"^ The 
description of this phenomena is performed by means of the fundamental magnitudes, 
energy, electron density and their derivatives with respect of the number of particles M 
which incorporates this change at the very basis of the descriptor deflnitionsd’^’^ It im¬ 
poses the knowledge of the energy 8^ dependence with this number, where M = N rkv 
with N E N and z/ G M the electron transferred fraction which lies in the u G (0, 1) 
interval (N and R fields of positive integer and real numbers, respectively), i.e., their 
dependence not only at integer numbers but for all values of Nevertheless the 

concept of reactivity is related to the interaction between the entities involved in the 
rearrangements, the quantities used to describe the behavior of reactive phenomena are 
usually evaluated by finite differences with respect to integer number of the particles 
of isolated speciesd’^ This approach ignores the values for the magnitudes at fractional 
numbers and consequently the electron exchange as the onset of chemical behavior. To 
avoid such inconsistences, exhaustive works supported by the above mentioned energy 
dependence has been performed about essential treatment for these chemical descrip¬ 
tors from state function approach.^ The most general way to describe a quantum state 
of a system is by means of the density matrix (DM) formalism®”® which contains 
the complete information about the system, i.e., from which all properties may be 
determined.®”® In the case of a system with a fractional number of particles, the state 
can not be described by a pure state neither a canonical ensemble but by an statistical 
ensemble of pure states with different number of electrons, i.e., a grand-canonical like 
ensemble (GC).®’^ For ground states the dependence of the energy 8^ and the DMs 
is a piecewise-continuous linear functions of A/” and only the bordering integers N and 
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iV ± 1 enter in this ensemble.^ Consequently, ground state properties then have the 
same dependence (cf. Fig. la of Ref. 3) and the hrst derivatives of the energy and the 
density are staircase functions of A/”, undehned at the integers and constant in between 
(cf. Fig. lb of Ref. 3). Thus, second derivatives vanish in between and are not dehned 
for integers.^ This behavior has deep consequences, for instance, electronegativity equa¬ 
lization principle does not apply and the electronic principles based on hardness has 
no rigorous foundations.^"^ As mentioned above, attempts to avoid the inconsistencies 
has been proposed for pure state formulations^ but no a general statistical formulation 
has been reported considering the system interaction with the environment, i.e., other 
subsystem or reservoir (S-R) interactions. The main objective of this report is to intro¬ 
duce the rigorous scheme for the DMs structure^’^ within the interaction between the 
subsystems to highlight the essence of the reactivity descriptors. 

The weighted sum of the complete set of all accesible M-electron pure state 
density matrices®"® ~ 1*^^ X *^^1 mixture, with |<F^ > the k — 

th quantum state function in the antisymmetric M-electron Hilbert space Jbu, 
represents the state of the system by the density matrix D, whose carrier space is the 
Fock space = 0^=0*^^ where 0 symbol indicates direct sum^’®, and reads 

11=^^ u^M |<Ff >< I; X] X] = 1; > 0 (1) 

M>0 fc>0 M>0 k>0 

where M > 0] k > 0} stands for the set of statistical weights, i.e., probability of 

occurrence of the pure state |<F^ > in the mixture. This state admits particle number 
fluctuation and hence the system may posses a non-integer number of particles. We will 
refer it as the GC distribution. D is an Hermitian, positive semi-dehnite (its eigenvalues 
are nonnegative or null), bounded (its elements are bounded) and hnite trace (sum of 
diagonal elements) matrix. Thus, because of its probabilistic interpretation it may be 
normalized to unity, i.e., tr{D) = 1.®"® 

The fundamental chemical concepts are the summary of the physical information 
contained in D and described by two types of quantities called descriptors. On one 
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side are those coming from the direct integration of a function of the density related 
to the classical concepts of chemistry® and on the other side, those related to the 
energy/density derivatives of first and higher ordersd Let us remark at this point 
that the physical domains within the molecular structure which possess a fractional 
number of particles M = N described by GC states (Eq. (1)) may be interpreted 
as an average.The convex structure of D for ground states may be expressed in 
two branches each one as a two state model of N and ± 1 Hilbert spaces as D = 
(1 - p) ^Dq + p with ^Do = 1$^ >< 4>^| and >< the 

corresponding ground state DMs. The energy is = {1 — p) Sq + p The 

onset of the above mentioned inconsistencies are related to the lack of information in 
the description about the interaction of the system (subsystem) with the environment. 
To avoid it we introduce a driven interaction potential Ui, which has its origin in the 
subsystem fragment within the Atoms in Molecules (AIM) framework^ or reservoir 
interactions effects® describing the influence of the environment. The subscript in JJ^, 
indicates its electron fraction dependence. So that, the Hamiltonian for each of the two 
branches, H = Hq + in matrix form is "H = >< >< 

$^^^1 + >< >< ^^1 where Hq stands for the isolated 

system Hamiltonian whose states produce discontinuities in the high order descriptors.^ 
The interaction potential induces a new solution D which at equilibrium exchanges 
electrons at a time independent rate, i.e., the electron exchange is constant in time, 
and hence the solution may exhibit a coherent structure® expressed as 

D = D + >< + A^* >< 4>^| (2) 

where the hrst term of the r.h.s. represents the solution for while the last two 

terms correspond to the coupling interaction of the > and > states. Within 

this scenario, the energy of the system under the influence of the environment results 

= Tr{nD) = + 2Me(W^^*A^) 


(3) 
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where the symbol Me stands for the real part of the complex number Therefore, 

this term may introduce a //-nonlinearity dependence for the energy and the DM and 
thus enable us to perform the calculation of the chemical descriptor of arbitrary order 
avoiding the discontinuity problem. To be more concrete let us write two important 
descriptors to show these ideas. The chemical potential dehned by^“"^ 





(4) 


at constant v (external held), becomes 


= 


2 


/Me(W±*A±) 




dv 


(5) 


where stands for the = —EA and = —IP, electron affinity and ionization 
potential, respectively and the second term in Eq. (5) shows the interaction contribution 
to the chemical potential which avoids the discontinuity^ and admits the equalization 
principle to be fulhlled^. For that goal, let us consider two fragments and VLb within 
a molecular framework which at equilibrium must obey the condition i-e., 

the chemical potential of the donor fragment must be equal to that of the acceptor 
fragment and it is the second term of the r.h.s. of Eq. (5) which enable this condition. 
The hardness^ which would vanish because of the chemical potential discontinuity for 
an isolated system, i.e., without interaction with an environment^, becomes non null 
due to the S-R interaction and reads 


1 (^ , ( d^ne(uri^i) \ 

^ 2 (aw j av ) 


( 6 ) 


Let us hnally to mention some concluding remarks. This approach is general regarding 
the DMs may be calculated from any arbitrary methodology and it depends only of the 
model used for describing the interaction with the environment; also this formulation 
permits to recover the piece-wise dependence as the interaction vanishes, i.e., \Jv —> 0 . 
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